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Q ' Abstract 

^ : 

The neural network approach is applied to search for the Z'-boson 
\ within the LEP2 data set for e^e~ — > e^e~ scattering process. In 

^> ' the course of the analysis, the data set is reduced by 20 percent. The 

■ axial-vector Cg and vector v\ couplings of the Z' are estimated at 
\ 95% CL within a two-parameter fit. The mass is determined to be 

0.53 < rrizi < 1-05 TeV. Comparisons with other results are given. 

(N 
O . 

g ; 1 Introduction 

A lot of extended models includes a so-called Z' gauge boson - massive neutral 
^ ' vector particle associated with an extra subgroup of underlying group 

■ [2]. In particular, this particle is predicted by numerous grand unification 
models. In all the models, the Abelian Z'-boson is described by a low- 
energy f/(l) gauge subgroup originated in some symmetry breaking pattern. 
Searching for Z' is one of the goals of future experiments at LHC [4J, and 
current ones at Tevatron. It can manifest itself either as a real or a virtual 
state dependently on the value of its mass. 

With this goals keeping in mind, it is reasonable to take into consideration 
all accessible nowadays information about Z'-boson, following from different 
experiments established at low energies. In particular, in LEP searching for 
Z' it was established a model dependent approach and low bounds on its 
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mass mz' > 400 GeV and/or m^/ > 800 GeV have been obtained [H], which 
are dependent on a specific model. No actual signals of the Z' have been 
detected. In recent experiments at Tevatron the derived low bounds are a 
little bit larger and correspond to the masses mz' > 850 GeV |4j. 



On the other hand, in series of papers [T] [2] a model independent method 
of searching for the Z' was developed. This approach accounts for a renor- 
malizability of an underlying unknown in other respect theory. The require- 
ment of renormalizability results in a series of relations between low energy 
parameters describing interactions of the Z' with fermions of the standard 
model. This reduces the number of low-energy parameters which must be 
fitted in experiments. In this way the couplings of Z' to the fermions of the 
standard model and the mass m^/ were estimated at 1 — 2 cr CL in the one- 
parameter [T] and many-parameter |2] fits. If these relations are not taken 
into consideration, no signals (hints, in fact) follow. 

It was also concluded that the statistics of the LEP experiments is not 
too large to detect the Z' as the virtual state with enough high precision. So, 
some further analysis is reasonable. One needs in the estimates which could 
be used in future experiments. To determine them in a maximally full way 
we address to the analysis based on the predictions of the neural networks 
(for applications in high energy physics see, for example, [7]). The main 
idea of this approach is to constrain a given data set in such a way that an 
amount of the data is considered as an inessential background and omitted. 
The remaining data are expected to give a more precise fit of the parameters 
of interest. This is the goal of the present investigation. We treat the full 
data set on the Bhabha scattering process obtained in LEP2 experiments by 
using the neural network method in order to determine the couplings to the 
fermions and the mass of the Z' boson. 

2 Neural network (NN) predictions 

The lack of statistics in LEP experiments does not allow the determination 
of the Z'-boson mass with CL more than 2a in one-parameter fit [T] and 
more than la in two-parameter fit [2]. We propose two points to overcome 
this restriction for the case of the Bhabha scattering process when a many 
parameter fit is applied. 

First, an increase in parametric space could be compensated by an in- 
crease in the data set, if all possible information is included in considera- 



2 



tion. In this paper we take into consideration the differential cross-sections 
measured by the L3 Collaboration at 183-189 GeV, OPAL Collaboration at 
130-207 GeV, ALEPH Collaboration at 130-183 GeV and the cross-sections 
obtained by the DELPHI Collaboration at energies 189-207 GeV. They form 
the set of differential cross-sections and their uncertainties for the center of 
mass energies. 

Second, we use NNs to predict the data set of the investigated process 
with increased statistics. 

In experiments, the Air scattering angle is divided into bins where the 
detectors are placed. For each bin we have a differential cross-section and 
an uncertainty [5]. The Z' extends the SM and therefore contributes to the 
differential cross-section 

da _ da^^ ^ da^ 
dz dz dz ' 

where z = cos9, 9 is a scattering angle, is the contribution coming due 

to the SM particles and is the contribution due to the Z' presence. In 
actual calculations, the first term was taken from the results reported by the 
LEP Collaborations and the second one has been calculated in the improved 
Born approximation with the relations due to renormalizability been taken 
into account. In this case it is possible to construct the observables which 
uniquely pick out the Z' virtual state [1],[2]. 

Hence, the signal of the Z' could be searched in the deviation of experi- 
mental differential cross-section from the SM differential cross-section. 

Taking into account all the noted above we prepare the needed exper- 
imental data set. It consists of a differential cross-section for the Bhabha 
process and an uncertainty. In order to use it in our analysis we subtract 
from it the SM differential cross-section for the Bhabha process. Going this 
way we obtain the data set: 

1. scattering sector (bin). 

2. - the difference between a measured differential cross-section and 

dz 

the corresponding SM differential cross-section. 

3. e - uncertainty. 

An overall distribution of the data is shown in FigH] As it is seen, this 
forms the set of segments determining the value of the cross-section with its 
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Figure 1: Experimental data are segments resulting from the experimental 
uncertainties. Bins are enumerated. 1-5 bins correspond to scattering in the 
backward direction {CosO < 0), the others - forward direction (Cos6 > 0). 

uncertainty. To obtain a more strict constraint on the axial-vector and 
vector v"^ couplings squared, we have to decrease the lengths of segments. 
For this purpose the NN analysis is applied. 

At first stage, the NNs were trained to recognize the signal and the back- 
ground. 

The contribution of the Z' to a differential cross-section for the Bhabha 
process obtained within the SM extended by Z' was taken as the signal. 
Analytic expression for this contribution reads [1]: 



where Fa , Fy,Fav are the functions depend on the SM couplings. In actual 
analysis they have been calculated with accounting for the relations due to 
renormalizability. This contribution was computed for the values of 0.0 < 
■y^ < 4 X 10~^ and 0.0 < < 4 x 10~^. Such the choice is based on the 
results, obtained in [Ij. 

As the background we set the deviation from the signal equaled to the 
redoubled uncertainty of LEP2 experimental data. Thus, the network was 
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trained to discard data, which correspond to large deviations from the signal, 
but include the ones corresponding to the probable signal of the Z' . 

The NN processing makes the cutoff of the data. The cutoff eliminates 
the data which are assumed to be the background signal. Side by side with 
the background signals the NN could discard the signal data too. But total 
amount of these events is negligibly small as compared to the amount of 
discarded background points. 

To create and train NNs the MLPFit program [6] was used. Three-layer 
NNs were used, with 2 neurons in input, 10 neurons in hidden and 1 neuron 
in output layer. 

An input vector for the networks consists of the scattering sector and the 
differential cross-section for this sector. The training algorithm with back 
propagation of errors was used. The type of training - with tutor - was 
applied. We also worked out a necessary computer program to solve the 
problem. 

During processing, the NNs discarded all the data that have produced a 
less than 0.9 at NN output. After NN processing, the length of segments 
was decreased. The total decrease of them is 3-27 percent (these values vary 
for different data sets). The comparative general plot of data set before and 
after the application of the NN is shown in Figl2l 

After NNs processing, the data were analyzed by means of the method. 
Denoting an observable by (jj, one can construct the x^-function 

X (a,^^e) = \_\ 7- ] , (3) 

i ^ 

where Aa*^^^ and 5cr is the deviation of the measured cross-section from 
the SM one and the uncertainty for the Bhabha process, and Ao"*'^ is the 
calculated Z' contribution. The sum in Eq.(l3]) refers to either the data in one 
specific process or the combined data for several processes. By minimizing 
the x^-function, the maximal-likelihood estimate for the Z' couplings was 
derived. The confidence area in the parameter space (w, a) corresponding to 
the probability /5 is defined as 

< xLn + XcL,/3 (4) 

where Xcl,/3 /3-level of the x^-distribution with 2 degrees of freedom. 

The analysis of the x^-function for the xP' observable gives Xmm — 158.49 
at = 2.37 X 10"^. The 95% CL area (x^l = 5.99) is shown in FigJSl From 



5 



0.2 



-0.2 



Bin 



10' 



11 



12 



Initial data 
Processed data 



Figure 2: Plot of e^e e+e process data before and after NN processing 



this area the value v'^ = (2.4 ± 1.99) x 10~^ is obtained. The 95% CL area 
for data before the NN processing is shown in FigjH 

The analysis of the x^-function for the observable results in Xmin — 
629.66 at = 1.6 x 10"^. The 95% CL area (xcl = 5.99) is shown in FigJSl 
To compare the results, the 95% CL area for data before NN processing is 
shown in FigJHl As one can see, the zero point of d is inside the confidence 
area. So, we are able to obtain only an upper limit on its value < 1.1 x 10~^ 
at the 2a CL. 

To constrain the value of Z' mass by the derived bounds on the four- 
fermion couplings v"^ let us assume that the coupling g is of the order of the 
SM gauge couplings, g'^/{47i) ~ 0.01 — 0.03. Then the obtained value 
corresponds to rriz' = 0.53 — 1.05 TeV. 

In Ref. [1] the value v'^ = 2.18 ± 1.82 x 10""^ was obtained for the Bhabha 
process at the 2o" CL within one-parameter fit. 



3 Discussion 

It is shown in Figl2]how the NN selects data corresponding to the background. 
We note once again that the NN discards not only a background but also 
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Table 1: Constraints on the Z' four-fermion couplings obtained in different 
investigations 



one-parameter fit 


two-parameter fit 


NN two-parameter fit 


1(7 


2a 


2a 


2a 


xlO-^ 


2.24 ±0.92 


2.18 ± 1.82 


< 1.44 


2.4 ±1.99 



the signal. Nevertheless, the discarded signal data are negligibly small as 
compared to the discarded background ones. 

The obtained values of = (2.4 ± 1.99) x 10"^ and < i.i x 10"^ 
are comparable to the results obtained in [Ij. In this paper the results were 
obtained within the one-parameter fit. In paper |2] it was shown that there 
is only a la CL hint in the two-parameter fit. The analysis carried out in the 
present investigation shows the 2a CL hint of the Z'. The estimate of the 
mass value gives mz' = 0.53 — 1.05 TeV. This is signaling a comparably not 
heavy Z', that is of interest for the experiments at the Tevatron and LHC. 
The estimated couplings are also important to analyze the current and future 
experimental data. 

To conclude, we note that the derived predictions of the NNs analysis 
of the two parameter fits are in good agreement with other one-parameter 
model- independent fits accounting for the Z' gauge boson existence [5J, [2]. 
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Figure 4: The 95% CL in v — a plane for f ^ observable. The data set before 
the NN processing. 
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Figure 5: The 95% CL in v — a plane for observable. The data set after 
the NN processing. As it is seen, the zero point of a is inside the confidence 
area. 
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Figure 6: The 95% CL m v — a plane for a? observable. The data set before 
the NN processing. 
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